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We introduce a spin ladder with antiferromagnetic Ising ZZ interactions along the legs, and inter- 
actions on the rungs which interpolate between the Ising ladder and the quantum compass ladder. 
We show that the entire energy spectrum of the ladder may be determined exactly for finite number 
of spins 2N by mapping to the quantum Ising chain and using Jordan- Wigner transformation in in- 
variant subspaces. We also demonstrate that subspaces with spin defects lead to excited states using 
finite size scaling, and the ground state corresponds to the quantum Ising model without defects. 
At the quantum phase transition to maximally frustrated interactions of the compass ladder, the ZZ 
spin correlation function on the rungs collapses to zero and the ground state degeneracy increases by 
2. We formulate a systematic method to calculate the partition function for a mesoscopic system, 
and employ it to demonstrate that fragmentation of the compass ladder by kink defects increases 
with increasing temperature. The obtained heat capacity of a large compass ladder consisting of 
2N = 104 spins reveals two relevant energy scales and has a broad maximum due to dense energy 
spectrum. The present exact results elucidate the nature of the quantum phase transition from 
ordered to disordered ground state found in the compass model in two dimensions. 
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I. INTRODUCTION 

Spin ladders play an important role in quantum mag- 
netism. Interest in them is motivated by their numer- 
ous experimental realizations in transition metal oxides^ 
and has increased over the last two decades. One 
of recently investigated realizations of spin ladders are 
Sr„_iCu„+i02n cuprates (with n = 3, 5, 7, • • and the 
simplest of them, a spin ladder with two legs connected 
by rungs, is realized in Sr2Cu406. Excitation spectra of 
such antiferromagnetic (AF) spin ladders are rich and 
were understood only in the last decade. They con- 
sist of triplet excitations, bound states and two-particle 
continuum^i and were calculated in unprecedented detail 
for quantum AF spin 5 = 1/2 two-leg ladder employ- 
ing optimally chosen unitary transformationii In some of 
spin ladder systems charge degrees of freedom also play a 
role, as for instance in Q:'-NaV2 05, where AF order and 
charge order coexist in spin ladders with two legSf^ or in 
the Cu-0 planes of Laa;Sri4_a;Cu2404i, where spin and 
charge order coexist for some values oi This advance 
in the theoretical understanding of the ground states and 
excitation spectra of spin ladders is accompanied by re- 
cent experimental investigations of triplon spectra by in- 
elastic neutron scattering^ of almost perfect spin ladders 
in La4SrioCu2404i. Finally, in the theory spin ladders 
could serve as a testing ground for new (ordered or dis- 
ordered) phases which might arise for various frustrated 
exchange interactionsi^ 

A particularly interesting situation arises when frus- 
tration of spin interactions may be tuned by varying 



strength of certain coupling constants, and could thus ex- 
hibit transitions between ordered and disordered phases. 
On the one hand, periodically distributed frustrated Ising 
interactions do not suffice to destroy magnetic long-range 
order in a two-dimensional (2D) system, but only re- 
duce the temperature of the magnetic phase transitioni^ 
On the other hand, when the model is quantum, in- 
creasing frustration of exchange interactions may trig- 
ger a quantum phase transition (QPT), as for instance 
in the one-dimensional (ID) compass modeliiS Physical 
realizations of frustrated interactions occur in 2D and 
three-dimensional spin-orbital models derived for Mott 
insulators in transition metal oxides in the orbital part 
of the superexchange. In such models frustration is in- 
trinsic and follows from the directional nature of orbital 
interactionsFii Usually such frustration is removed ei- 
ther by Hund's exchange Jh or by Jahn- Teller orbital 
interactions, but when these terms are absent it leads 
to a disordered orbital liquid ground state. Perhaps the 
simplest realistic example of this behavior is the (Kugel- 
Khomskii) model for Cu^+ ions in (P electronic configu- 
ration at Jh = 0, where a disordered ground state was 
foundji^ Examples of such disordered states are either 
various valence-bond phases with singlet spin configura- 
tions on selected bonds^^ or orbital liquids established 
both in tig systemsi^ and in Cg systemsi^^ Characteris- 
tic features of spin-orbital models are enhanced quantum 
effects and entanglement so their ground states can- 
not be predicted using mean-field decoupling schemes. 
Also in doped systems some unexpected features emerge 
for frustrated orbital superexchange interactions, and the 
quasiparticle states are qualitatively different from those 
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arising in the spin t-J modelJi Therefore, it is of great 
interest to investigate spin models with frustrated inter- 
actions which stand for the orbital part of the superex- 
change, particularly when such models could be solved 
exactly. 

Although the orbital superexchange interactions arc 
frequently Ising-like, they lead to quantum models with 
intrinsically frustrated exchange models as different or- 
bital components interact depending on the bond orien- 
tation in real spacC]^ A generic case of such frustrated in- 
teractions is the so-called 2D quantum compass modclji^ 
which was recently investigated numericallyj ^ ^ Al- 
though orbital superexchange interactions in Mott in- 
sulators are typically AF j^^i^^i^^'^"^ a similar frustration 
concerns also ferromagnetic (FM) interactions, and a 
QPT was also found in the compass model with FM 
interactions 

The ID variant of the compass model with alternat- 
ing interactions of z-th and x-th spin components on 
even and odd bonds was solved exactly by an analyti- 
cal methodfiSi and entanglement in the ground state was 
analyzed recently]^ We note that the ID compass model 
(the model of Ref. [13 in the limit of equal and alter- 
nating interactions on the bonds) is equivalent to the ID 
anisotropic XY model, solved in the seventies An exact 
solution of the ID compass model demonstrates that cer- 
tain nearest-neighbor spin correlation functions change 
discontinuously at the point of a QPT when both types of 
interactions have the same strength. This somewhat ex- 
otic behavior follows because the QPT occurs at the mul- 
ticritical point in the parameter space A similar dis- 
continuous behavior of nearest-neighbor spin correlations 
was also found numerically for the 2D compass modcljii^ 
While small anisotropy of interactions leads to particular 
short-range correlations dictated by the stronger interac- 
tion, in both ID and 2D compass model one finds a QPT 
to a highly degenerate disordered ground state when the 
interactions are balanced. 

The purpose of this paper is to present an exact solu- 
tion of the compass model on a spin ladder, with ZZ Ising 
interactions between z-th spin components along the lad- 
der legs, and interactions on the rungs which gradually 
evolve from ZZ Ising interactions to XX Ising ones. In 
this way the interactions interpolate between the clas- 
sical Ising spin ladder and the quantum compass ladder 
with frustrated interactions. The latter case will be called 
compass ladder below — it stands for a generic competi- 
tion between orbital interactions on different bonds and 
can serve to understand better the physical consequences 
of the frustrated orbital superexchange. 

The paper is organized as follows. The model and 
its invariant dimer subspaces are introduced in Sec. II. 
Next the ground state and the lowest excited states of 
the model are found in Sec. Ill by solving the model in 
all nonequivalent subspaces. Thereby we discuss the role 
played by defects in spin configuration and show that the 
ground state is obtained by solving the ID quantum Ising 
(pseudospin) model (QIM). Using an example of a finite 



system, we provide an example of the energy spectrum, 
and next extrapolate the ground state energy obtained 
for finite systems to the thermodynamic limit. We also 
present the changes of spin correlations at the QPT, and 
derive the long-range spin correlations. Next we con- 
struct canonical ensemble for the spin ladder in Sec. IV 
and present the details concerning the calculation of ener- 
gies in the appendix. The constructed partition function 
is used to derive such thermodynamic properties of the 
compass ladder as the temperature variation of spin cor- 
relations, and the average length of fragmented chains 
separated by kinked areas in Sec. V. In Sec. |VT] we 
present the evolution of heat capacity Cy when inter- 
actions change from the Ising to compass ladder for a 
small ladder of TV = 8 spins, and next analyze Cy for 
a large (mcsoscopic) compass ladder of 2N = 104 spins. 
While the characteristic excitation energies responsible 
for the maxima in heat capacities can be deduced from 
the energy spectrum for N = 8 spins, generic features 
of excitations follow from the form of Cy in case of the 
mcsoscopic compass ladder. Final discussion and sum- 
mary of the results are given in Sec. VII. 

II. COMPASS MODEL ON A LADDER 

We consider a spin ladder with N rungs {2i — l,2i) 
labelled by i ~ 1,2, ■ ■ ■ , N . The interactions along lad- 
der legs are Ising-like with AF coupling J between z- 
th spin components (cf (7f_|_]^), while AF interactions on 
the rungs interpolate between the Ising coupling of z-th 
(2crfj_]^o'^_l_]^) and x-th {2a^_i(j'^^i) spin components, 

N 

n{a) = 2j5^{«a2Vi4, + (l-«)'^2Vi'^2j 

i=l 
N 

+ (a2^_i(Tf,+i + «+2) , (2.1) 

i=l 

by varying parameter < a < 1. We assume pe- 
riodic boundary conditions along the ladder legs, i.e., 

(y2N+i = '^1 ^^"^ ^2Af+i = '^2- The factor of two for the 
interactions on the rungs cx 2 J was chosen to guarantee 
the same strength of interactions on the rungs (with only 
one rung neighbor of each spin) as along the ladder legs 
(with two leg neighbors). Increasing a gradually modifies 
the interactions on the rungs and increases frustration. 
For a = one finds the reference Ising ladder, while at 
a = 1 the interactions describe a competition between 
frustrated ZZ interactions along the ladder legs and 2XX 
interactions on the rungs, characteristic of the compass 
ladder. A representative compass ladder with N ~ 4 
rungs (i.e., 2N = 8 spins) is shown in Fig. [TJ 

To solve the spin ladder given by Eq. (|2.ip in the 
range of < a < 1 we notice that [H(a), (T|j_]^(T|j] = 0. 
Therefore we have a set of N symmetry operators, 

R, = ol,_^ol, , (2.2) 
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with a constant 



FIG. 1: (Color online) Schematic view of the quantum com- 
pass ladder with A'' = 4 rungs, described by Hamiltonian (|2.ip 
with a — 1. Interactions along the ladder legs labeled as 
ZZ (horizontal lines) are (^2i-i'^2i+i (upper leg) and o"2if"2i+2 
(lower leg) . The interactions along the rungs labeled as 2XX 
(vertical lines) are 2crfj_i(jfj (the factor of 2 simulates the 
periodic boundary condition along the rungs). Dashed lines 
indicate periodic boundary conditions along the ladder legs. 



with respective eigenvalues ~ ±1. Each state of the 
system can be thus written in a basis of af eigenvectors 
|si7 ^2, S3, . . . , S2n) fixed by strings of quantum numbers 
Si = ±1. These vectors can be parametrized differently 
by a new set of quantum numbers {i,} and {ri}, with 
i = 1, 2, • • • , A^; they are related to the old ones by the 
formulae: ti = S2i-i and = S2i-iS2i- Now we intro- 
duce new notation for the basis states 



1*1, ^2, ■ • ■ , *Af)rir2---r„ = 1*1 > ^l''! : *2 , *2''2, • ■ • , *Ar, ^WAf) , 

(2.3) 

where the right-hand side of Eq. (|2.3p is the state 
S2, S3, . . . , S2Ar) Written in terms of variables {ti} and 
{vi}, and the left-hand side defines new notation. This 
notation highlights the different role played by r^'s, which 
are conserved quantities, and by i^'s, being new pseu- 
dospin variables. For states like in Eq. (|2.3p . we define 
new pseudospin operators and acting on {ti} quan- 
tum numbers as Pauli matrices, e.g. for i = 1: 



Tl\tl,t2, ■ ■ ■ , *Af) 

Tl 1*1, *2, ■ ■ ■ , tN)rir2---rN 



I ~ *1, *2, ■ • ■ , tN)rir: 



ti 1*1, *2, . . . , *Ar)j-i 



T-2---rjv ■ 



(2.4) 



A similar transformation was introduced for a frustrated 
spin-1/2 chain by Emery and Noguera^^ who showed 
that it can by mapped onto an Ising model in a trans- 
verse field. Recently this procedure was used to investi- 
gate quantum criticality in a two-leg strongly correlated 
ladder model at quarter fiUingi^l 

The Hamiltonian can be now written in a common 
eigenbasis of Ri (|2.2p operators by means of {rf , r/} op- 
erators. In a subspacc labelled by a string ri, r2, ■ ■ ■ , r^v, 
the reduced form of the Hamiltonian is 



Tin 



N 

r^.-.r^ia) = J^{(I + r,r,+i)rfTf+i +2aTf} 
1=1 

+ 2Ja-(") , (2.5) 



N 



Cp{a) = (1 - g) 



N+l 



(2.6) 



This leads 



and periodic boundary condition r 
to the exactly solvable QIM with transverse field )^^'^^i'^° 
if only ri = 1 or = —1. Otherwise there are always 
some t^tI^i interactions missing (defects created in the 
chain) and we obtain a set of disconnected quantum Ising 
chains with loose ends and different lengths. The bonds 
with no pseudospin interactions may stand next to each 
other, so in an extreme case when fi+i = — for all i, 
one finds no Ising bonds and no chains appear. 

One may easily recognize that the ground state of the 
spin ladder described by Hamiltonian (|2.ip lies in a sub- 
space with r,; = —I for a < I. First of all, = —I 
minimizes Cp{a), see Eq. (|2.6p . To understand a second 
reason which justifies the above statement let us examine 
a partial Hamiltonian (open chain) of the form 



L-l 



L 

t^,+2JaJ2rr, (2.7) 

4=1 



with 2 < i < — I. Note that it appears gcnerically 
in Eq. (j2.5p and consists of two terms containing pseu- 
dospin operators {rf } and {rf }. Let us call them 
and and denote the ground state of as \x) with 
energy iS^,. The mean value of 7i(a, L) in state |.t) is also 
Ex because every rf operator has zero expectation value 
in state jx), i.e., (xjrf |x) ~ 0. However, wc know that 
|a;) is not an eigenvector of 7i(a, L) which implies that 
7i(a, L) must have a lower energy than in the ground 
state. This shows that the presence of rf rf^j^ bonds in the 
Hamiltonian 7i(a, L) lowers the energy of bare H^. One 
may also expect that this energy decreases with increas- 
ing length L of the chain, and is proportional to L in the 
thermodynamic limit. The numerical evidence for this 
are plots of the ground state energy versus L presented 
in section 3. Looking at Hamiltonian (|2.5p we see that 
the longest chains of the type (|2.7[) appear in subspaces 
with ri = —1 and = I, but the constant term Cf(a) 
favors = — 1 if only a < 1. For a = 1 the ground state 
can be in both subspaces, and its degeneracy follows, see 
below. 



III. ENERGY SPECTRA IN INVARIANT 
SUBSPACES 

A. Quantum Ising model 

To find the ground state of spin ladder (|2.ip wc need to 
solve the QIM that arises from Eq. (|2.5p when = — 1. 
Thus we need to diagonalize the Hamiltonian of the form 



N 



Hqim(/3,«) = 2J^(/3r;T4i +aTf) , 



(3.1) 
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which is related to our original problem by the formula 



n-1-i...-i = 7^qim(1, a) - 2NJ{1 - a) 



(3.2) 



The formal parameter /3 is introduced for convenience 
and will be used to determine the correlation functions 
along the ladder legs by differentiation, see below. The 
standard way of solving TYqim starts with Jordan- Wigner 
(JW) transformation. This non-linear mapping replac- 
ing spin operators by spinless fermions is of the form 



(1 



(3.3) 



The boundary condition for fermion operators {q} after 
inserting them into TYqim (|3.ip is antiperiodic for even 
and periodic for odd number of JW quasiparticles in the 
chain. The operator V of the parity of fermions, 



N 



i=l 



(3.4) 



corresponds to the operation of flipping all spins along 
the z-th axis and commutes with TYqim- Therefore, the 
Hamiltonian can be split into two diagonal blocks, for 
even (+) and odd (— ) number of JW fermions by means 
of projection operators ^(1 ± V). Therefore we write 



1 



1 



(3.5) 



where 



N 

E 

1=1 

2JNa 



- Ci){cl^^ + c,+i) - 2ac.|c,;| 



(3.6) 



with two different boundary conditions: Ctv+i = Tci for 
(±) subspaces. Let us point out that the only conse- 
quence of the nonlinearity of the JW transformation is 
the minus sign which appears in the first bracket multi- 
plying p. This is thanks to one-dimensionality and only 
nearest-neighbor interactions in the reduced Hamiltonian 
(|2.5p . but is not the case for the original Hamiltonian 

(EH). 

Next step is the Fourier transformation. 



1 



(3.7) 



with quasimomenta k = ±(2/ - l)7r/iV = 1, 2, • • • , N/2] 
in an even subspace (-I-), and k ~ Q,7r,zL2lTr/N [I ~ 
1, 2, • • • , (A^/2 — 1)] in an odd one (— ). After transform- 
ing the operators in Eq. p.6p we obtain in a block 
diagonal form, 



4J 



2J 



k 

— ^- 

k 



{(3 COS k — a)c^Cfc 



+ h.c.) + 2JNa. (3.8) 



Diagonalization is completed by a Bogoliubov transfor- 
mation, defining new fermion operators 7^ = afcc[. -I- 
PkC-k (for A; ^ 0, TT, while the operators cq and c-^ have 
no partner and are left untransformed). Transformation 
coefficients ak and Pk arc obtained from the condition 



'Hqim, -fk 



uJkll 



(3.9) 



which is an cigenproblem in linear space spanned by op- 
erators and c_fe. We get two eigenvectors {ak,f3k), 
corresponding to the quasiparticle operators 7^ and 7-fc, 
and two corresponding eigenvalues u>k = ^Ek, with 



Ekif3,a) = Aj{a'^+ 13"^ ~2a(3 cos k} 



1/2 



(3.10) 



Therefore, the Hamiltonian is brought to the diagonal 
form in both subspaces 

n+ = E^^'^- - 1) , (3.11) 

n- = Y.'Ek (^jhk - + 4J(/3 - a)4co 

- 4J(/3 + a)4c^ + 4Ja . (3.12) 

We still need to transform the parity operator T'. Luck- 
ily, the Fourier transformation does not change its form 
and to see that so does the Bogoliubov transformation, 
one can look at the vacuum state |0) for quasiparticle 
operators 7^,. From the condition 7fc|0) = for all k we 
get 



10) = n 



ak + PkC_kc\ 



(3.13) 



where \vac) is a true vacuum state for JW fermions or 
a state with all spins up. From the form of |0) we see 
that it contains a superposition of all even numbers of 
quasiparticles cj^, and the total quasiparticle number is 
not fixed. Acting on the vacuum with a single creation 
operator 7]^ we obtain a state with odd number of JW 
fermions, because 7^ is a linear combination of a creation 

and annihilation c-k operator of a single fermion. In 
this way one may get convinced that the parity of quasi- 
particles 7jj! and is the same. 



B. Ground state and the energy spectrum 

From the diagonal form of the QIM Hamiltonian given 
by Eq. (|3.1ip we see that the ground state of spin ladder 
(|2.ip is simply |0) in subspace = — 1 (or = 1 when 
a = 0). For the ground state energy, one uses Eq. 
to get 



£^QiM(l,a)-2iVJ(l-a) 



(3.14) 
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FIG. 2: Eigenenergies En of the spin ladder (|2.ip of Fig. [T] 
with = 4 rungs for increasing q, obtained by exact diag- 
onalization. Different panels show energies in invariant sub- 
spaces of the effective Hamiltonian (|2.5p . with 1 and 1 stand- 
ing for positive or negative values of r;: (a) 1111, (b) llll, 
(c) 1111, (d) 1111, (e) 1111 and (f) 1111. While the sub- 
spaces (a) and (f) are unique, other subspaces are equivalent 
by symmetry to those shown in panels (b)-(e), resulting in 
total spectrum of 256 eigenstates. Quantum phase transition 
occurs at a = 1, where the lowest eigenenergies in the sub- 
spaces (a) and (f) become degenerate. In the thermodynamic 
limit N ~* oo the spectrum changes qualitatively — the two 
lowest energies in the subspaces 1111 and 1111 are degenerate 
and the ground state from the subspace (b) (111!) becomes 
the first excited state of the spin ladder. 



with -Eqim(1, ct) given in the thermodynamic limit by an 
integral 



£^qim(/?, a) 



N 
'2^ 



dkEk{l3,a). 



(3.15) 



The ground state in the absence of transverse field (at 
a = 0) is doubly degenerate — it is given by two possible 
Neel states. At finite a > 0, this degeneracy is removed, 
and the sum of the two Neel states (symmetric state), 
|0-|-), is the ground state, while their difference (antisym- 
metric state) becomes the first excited state. This first 
excited state, |0_) = 7JIO+), stems from the same sub- 
space and belongs to the spectrum of 7i~. The splitting 




FIG. 3: (Color online) Nearest neighbor correlation functions 
in the ground state for spin ladder (|2.ip in the thermodynamic 
limit A'^ 00. For increasing a spin correlations (crfi-icfi) 
on the rungs decrease from zero to — 2/7r. At the same time, 
AF correlations (o"2i_iO"2i-i-r) along the ladder legs gradually 
weaken (increase from the classical value —1 at a = to —2/n 
at q: = 1), and become degenerate with the rung (crfi„icrfj) 
correlations at the quantum critical point a — 1. Correlation 
function (cr|i_i(T|j) on the rungs, directly related to the sub- 
space indices Vi, remains constant {{o'2i-iO'2i) ~ ~1) in the 
entire range of a < 1, and jumps to at q = 1. 



of the states |0_|.) and |0_) increases with a, see Fig. HJa). 
For finite N and a > there is always finite energy dif- 
ference between the energies of |0+) and |0_) = 7^10+) 
states. However, in the thermodynamic limit N —^ 00, 
this energy gap vanishes for a < 1 and starts to grow as 
A J a at a = 1. 

The full spectrum for the ladder with TV = 4 rungs be- 
longs to six classes of subspaces equivalent by symmetry 
— it is depicted in Fig. [51 With increasing a the spec- 
trum changes qualitatively from discrete energy levels of 
the classical Ising ladder at a = 0, with the ground state 
energy per spin equal —2 J, to a narrower and quasi- 
continuous spectrum when the quantum compass ladder 
at a = 1 is approached, with the ground state energy 
—AJ/tt per spin. At the a = 1 point one finds an ad- 
ditional symmetry; subspaces indexed by r and — r are 
then equivalent which makes each energy level at least 
doubly degenerate. 



C. Correlation functions 

All the nontrivial nearest neighbor spin correlation 
functions in the ground state can be determined by tak- 
ing derivatives of the ground state energy Eqim{P, a) 
(|3.15p with respect to a or /3, while the others are ev- 
ident from the construction of the subspaces. In this 
way one finds {o'2i-i'^2i+i) correlation along the legs and 
(^fi-i'^'fi) along the rungs, shown in Fig. [31 Spin cor- 



'2i-l"2i 

relations 



along the legs increase from the 
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classical value —1 up to —Ij-K for a ~ 1. By symme- 
try, both ladder legs are equivalent and {<J2i-i'^- 



tions (crfi_iCr' 



2i-l"2i+l/ ~ 

for a = a;, z. At the same time spin correla- 
2i) along the rungs gradually develop from 
in the classical limit to — 2/7r at the quantum critical 
point a = 1. Both functions meet at a = 1 which indi- 
cates balanced interactions — ZZ along the legs and 2XX 
along the rungs in case of the quantum compass ladder 
(see Fig. [1]). 
For the remaining correlations one finds 



= , 



(3.16) 
(3.17) 



Eq. (|3.16p follows from the fact that operators <72i-i'^2i+i 
do not commute with the symmetry operators Ri (|2.2p . 
In turn, averages of the symmetry operators along the 
rungs (|3.17[) are constant and equal —1 for a < 1, but at 
a = 1 they change in a discontinuous way and become 
zero, because at this point the degeneracy of the ground 
state increases to 2 x 2 = 4, and the spins on the rungs 
are disordered, so the ZZ correlations vanish. 

Finally, one can calculate the long range correlation 
functions for z-th spin components. 



(3.18) 



The right-hand side of Eq. (|3.18[) can be obtained from 
the QIM by the so-called Toeplitz determinant^ and 
can be also found in Ref. AH the long range XX 

correlation functions are zero in the ground state as they 
do not commute with i?i's operators (|2.2p . 

Note that correlations {t^t^) vanish in any subspace 
when \i— j\ exceeds the length of the longest Ising chain. 
This is due to the fact that, as already mentioned in 
section HIl the effective Hamiltonian in a given subspace 
describes a set of completely independent quantum Ising 
chains. Thus, at finite temperature, one can expect that 
the compass ladder will be more disordered than a stan- 
dard, ID QIM. The problem of chain partition at finite 
temperature will be discussed in detail below. 



D. Energies in the subspaces with open Ising 
chains 



As already mentioned, the general Hamiltonian of the 
form (|2.5p is exactly solvable only in cases when = r^+i 
or r,; = — r^+i for all i. Therefore, one may find exactly 
the ground state of spin ladder (|2.ip . see below). Oth- 
erwise, in a general case (i.e., in arbitrary subspace) one 
needs to deal with a problem of the QIM on an open 
chain of length L where L < N, described by Hamilto- 
nian (123; 



L-l 



(3.19) 




-2.6 



FIG. 4: (Color online) Ground state energies per site e(a, L) 
for the QIM on open chains (|2.7|l as functions of inverse 
cliain's length 1/L (for 2 < L < 61) plotted for a = 1, 3/4, 
1/2 and 1/4, from bottom to top. Linear fit gives the asymp- 
totic values of energies for L oo, indicated by crosses; these 
are: e(a,L -> oo) = -2.55, -2.29, -2.13 and -2.03 for the 
respective values of a. 



After applying the JW transformation (|3.3p . Eq. p.l9p 
takes the form 



nia,L) - 2JE{(c1-c,)(4i 



2JLa , 



'j+i) - 2aclci} 

(3.20) 



with an open boundary condition c]^_^_i = 0. This condi- 
tion prevents us from the plane waves expansion, but we 
can still use the Bogoliubov transformation. We remark 
that the broken chain considered here is sufficient to get 
a general solution, and the sum over all subspaces with 
open (broken) chains is included in the partition function 
Z{a), see Sec. WJ\ 



We define new fcrmion operators as follows 

L 



J2 ("y c] , 



(3.21) 



for i = 1,2, . . . , L. Coefficients and Pij can be chosen 
in such a way that the transformation is canonical and 
Hia, L) takes the diagonal form: 

n{a, L)^J2 ^^("' ^) (^h - I) ■ (3-22) 
i=i ^ ^ 

Both excitations energies Ei and transformation coeffi- 
cients {aij,Pij} can be determined from the condition 



[■Hia,L),^}] = Ea} 
This leads to an eigenequation 



A B 

-B -A 



(3.23) 



(3.24) 
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FIG. 5: (Color online) Coefficients a{a) (lower line) and b{a) 
(upper line) of the linear fit (|3.28|) performed for the points 
(1/60, e(a, 60)) and (1/61, e(a, 61)) for different values of a. 
At Q = one recovers the classical values of the Ising chain. 

where A and B are matrices of size L x L {A is a sym- 
metric and B is an antisymmetric matrix), and di, Pi 
are vectors of length L. The explicit form of A and B for 
L = 4 is 



2a 


1 








1 


-2a 


1 








1 


-2a 


1 








1 


-2a 



(3.25) 




FIG. 6: (Color online) Excitation energy A(J) as a function 
of Q for: (i) a ladder with A'' = 4 rungs (lower line) and 
(ii) an infinite ladder (upper line). In the first long 
as the QPT is not approached, the excited state remains in 
the ground state subspace with ri = —1 and contains one 
Bogoliubov quasiparticle with k — n. Excitation energy is 
small, starts from zero and decreases quickly with growing 
N. The latter excited state collapses to the ground state for 
A = oo, so the first excitation is here different than the one 
for infinite A. On the contrary, the excited state for N — oo 
contains two Bogoliubov quasiparticles with fc = 0*. This 
leads to the linear gap following A(q) = 8J(1 — a). 



and 



B^2J 



1 0\ 

-10 10 
0-101 
-1 o/ 



(3.26) 



which can be simply generalized to the case of any finite 
L. The spectrum of 7i(a,L) can be now determined by 
a numerical diagonalization of the 2L x 2L matrix from 
Eq. (|3.24p . For each L one obtains a set of 2L eigen- 
values symmetric around zero. Only the positive ones 
are the excitation energies Ei appearing in Eq. (|3.22p . 
Therefore, the ground state energy EQ{a,L) is obtained 
in absence of any excited states, so the energy per site 
can be easily expressed as 



1 1 ^ 

e{a,L) = -£'o(a,L) = - ^ ^E.i{a,L 



(3.27) 



Fixing a and increasing L wc can trace the depen- 
dence of e{a, L) on the system size and make an ex- 
trapolation to an infinite chain L — > oo. Results for 
e[a,L) (|3.27p as a function of decreasing 1/L, obtained 
for a = 1, 3/4, 1/2, 1/4 and L changing from 2 to 61, are 
shown in Fig. [H The energies decrease with increasing 
L which suggests that the ground state corresponds in- 
deed to a closed chain without any defects, as presented 
in Sec. [TnBl 



The dependence of e(a,L) on 1/L seems to be almost 
linear in each case. This is almost exact for a = 1 and for 
a — 1/4, while it holds approximately for intermediate 
values of a for in the regime of sufficiently large L. This 
observation can be used to derive a simple, approximate 
formula for the energy e(a, L). One can take the values 
of e(a, L) obtained for two largest L {L = 60, 61) with 
fixed a and perform a linear fit. Hence, we get 



e(a, L) = a{a) — — h{a) 
ij 



(3.28) 



with coefficients a and h depending on a. These new 
functions can be determined numerically for a chang- 
ing between and 1 with sufficiently small step. Re- 
sults obtained by a numerical analysis are plotted in Fig. 
m Both a and h starts from a value 2 at a = 0, then 
0(0;) decreases monotonically to about 0.72 while b{a) 
slightly increases to 2.55 ai a — 1. Eq. (|3.28p is exact 
for a = and any L, as well as for L — 60, 61 and any 
a. Nevertheless, looking at Fig. [H one can expect it 
to be a good approximation in case of sufficiently large 
L. From this formula one can read that for L ^ 00 one 
gets Eo{a,L) = —Lb{a) + 0(i°) which agrees with the 
classical intuition based on extensiveness of the internal 
energy. 
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E. Lowest energy excitations 



As we pointed out in See. IIII Bl the lowest exeited 
state in the case of a finite system, for a far enough from 
a = 1, is simply 7||0+) and belongs to the subspace 
Vi = —1. This is a collective excitation creating a wave 
of spin-flips in the ground state. Close to a = 1 one finds 
that the lowest excited state is the ground state from the 
subspace = 1 which means that the spin order along 
the rungs changes from AF to FM one along the z-th 
axis. 

The lowest energy excitation changes qualitatively in 
the thermodynamic limit N —> oo, where 7^10+) and 
|0+) states have the same energy and the dominating 
excitation is a pair of Bogoliubov quasiparticles with k ~ 
0^ which corresponds to fiipping one spin at a = 0. The 
first excited state remains in the = — 1 subspace for 
all a and the gap follows linear law A(a) = 8J(1 — a), 
see Fig. [6] This shows that in the thermodynamic limit 
{N oo) the low energy spectrum of the ladder is the 
same as for ordinary QIM. Note that such behavior is in 
sharp contrast with the case of finite ladder of = 4 
rungs. 



IV. CANONICAL ENSEMBLE FOR THE 
LADDER 



the lengths of these chains put in descending order. In 
case of our exemplary string its chain configuration is 
{4,4,3,3,2,2}. Variables {Li} must satisfy three condi- 
tions: (i) Li > 2 for all i, (ii) X^IIli — ^'^d (iii) 
n™i(-l)^' = (-!)"• The first two of them are obvi- 
ous, while the last one is a consequence of the periodic 
boundary conditions. Using chain parameters the effec- 
tive Hamiltonian TLrir2 - rN can be written as a sum of 
commuting operators 



m K 

^H{a,Li) - 2Ja Vrf 



i=l 



2JCp{a) 



(4.3) 

where K — N — stands for the total size of 

kinked areas. This formula refers to all subspaces ex- 
cluding those with r^-i-i = r^, where we have already ob- 
tained exact solutions. The evaluation of the constant 
Cp{a) can be completed by considering chain and kink 
areas in each subspace, see appendix. Having the diag- 
onal form of Ti.{a, L), given by Eq. (|3.22p . one can now 
calculate partition function for the ladder of 2A^ spins. 
It can be written as follows 

{Li}RiL,y 

+ Zo{a) , (4.4) 



A. Partition function 

In order to construct the partition function of spin lad- 
der (|2.1|) , we shall analyze its quantum states in different 
subspaces. Every invariant subspace introduced in Sec. 
im is labelled by a string rir2 ■ ■ ■ tn ■ Let us consider an 
exemplary string of the form 



1111111111111111111111 



(4.1) 



where 1 — —1, and either = TiJ^i or ^ ri+i. Each 
time when = r^^-i the chain continues, and when 
Ti ^ Ti-i-i we may say that a kink occurs at site i in 
the chain. We introduce a periodic boundary condition, 
so the string is closed to a loop and stands next to 
ri. From the point of view of the reduced Hamiltonian 
■■tn^ given by Eq. (|2.5p . it is useful to split the 
string {ri} into chains and kinked areas. A chain is a 
maximal sequence of ri 's without any kinks consisting at 
least of two sites. Kink areas are the intermediate areas 
separating neighboring chains. Using these definitions we 
can divide our exemplary string (|4.ip as follows 

1 1 1) 1 1 1 (1 1) (1 1 i) (1 1) i (1 1 1 1) (1 i 1) (1 , 

(4.2) 

where we adopt the convention to denote chains as 
{TiTi+i ■ ■ ■ Ti+p), and kink areas as )nri+i • • • ri+q{. For 
any string of ri 's containing m chains we can define chain 
configuration {Li} with i = 1, 2, . . . , m, where L^'s are 



where the sum over all {r} subspaces is replaced by 
sums over all chain configurations {Li{ and all R — 
'^i configurations possible for a given {Li{. Fac- 
tor i^Q[{Lj}, is a number of r subspaces for fixed 
chain configuration and fixed R when a < 1, and for 
a = 1 it is a number of r subspaces when only {Li{ 
is fixed. Partition function for any subspace containing 
open QIM chains or kinked areas is given by 



Z[a,{L,}) 



2^ cosh 



K 



2J 



Ta 



li 



2T 



(4.5) 



where {U} {i = 1,2, ... ,n) arc the different lengths of the 
chains appearing in the chain configuration {Li}, N(li) 
stands for the number of chains of the length li, and 
T is temperature in units of fc^ = 1. For example, the 
chain configuration {4, 4, 3, 3, 2, 2} of Eq. (|4.2p has n = 3, 
{k} {4,3,2} and N{k) = 2. The term Zo{a) is a 
contribution from subspaces with r^+i = r^. Using exact 
solutions p. lip , available in these subspaces, one finds 
that 



Zo{a) 



cosh 



2J 



N{1 - a) 



T 

/N-l N-1 

X E n ^ + ^ n ^i^i^ -f- '(4-6) 
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where the quasiparticle energies are: 

1/2 

E+ = 2j|l + a2 + 2acos (^^^^TT^l , (4.7) 

E- = 2j|l + a2 + 2acos ^^^^TT^l . (4.8) 

Appearance of both sine and cosine hyperboUc functions 
in Zq (|4.6p is due to the projection operators V intro- 
duced in section UlIAl 

B. Combinatorial factor 

To obtain numerical values of the partition function 
one has to get the explicit form of the combinatorial fac- 
tor Fa[{Li\,R^i^.T^]. This can be done in a simple way 
only for a = 1 when Cf(Q;) = 0, see Eq. (|2.6p . Then we 
have 

F^=i[{U},R{L^}]=Fi[{L,}] , (4.9) 

where Fi[{Li}] is the number of different r subspaces 
that can be obtained from a fixed chain configuration 
{Li\. Now we can derive a formula for this combinatorial 
factor. 

The chains can be put into the string in any order 
and these of equal length are indistinguishable. Apart 
from chains, there are also r^'s belonging to the kinked 
areas which determine the actual string configuration. 
We have K = N — X]t=i Li of them, they are indistin- 
guishable and can be distributed among m kinked areas. 
These degrees of freedom lead to a combinatorial factor 

N{h)\...N{l^)\\ K J' 

where li,l2, ■ ■ ■ ,ln {n < m) are the lengths of the chains 
without repetitions and N(li) is a number of chains of 
the length After determining the length of the first 
chain Li and the size of its kink area Ai, we still need 
to fix the position of ri. We have exactly Li + A\ possi- 
bilities. Next, we have to sum up over all possible values 
of L\ (which are Zi, ^2, • ■ • , ^n), all possible sizes of the 
kink area A\ (which are 1,2,..., and multiply by a 
combinatorial factor (|4.10p calculated for the remaining 
part of the string. The result is 

xf:«.+<.)(''--+_"r^).(4.ii) 

a=0 ^ ' 

where the factor of 2 in front comes from the fact that 
r\ = ±1. This number tells us how many times a given 
energy spectrum repeats itself among all subspaces when 
a = 1. The binomial factor appearing in formula (|4.1ip 
needs to be generalized with F functions when m = 1. 



V. COMPASS LADDER AT FINITE 
TEMPERATURE 

A. Correlation functions and chain fragmentation 

Nearest-neighbor correlation functions can be easily 
derived at finite temperature from the partition function 
Z(a, /3, 7), if we substitute our initial Hamiltonian 7Y(q:) 
given by Eq. (|2.ip by 

AT 

W(a, /?, 7) = 2J^ {l<yl^-xOl^ + (1 - 
i=l 
N 

+ J^Y. + «+2) ■ (5.1) 

i=l 

Then, after calculating the partition function, we recover 
spin correlations by differentiating Z{a,i3,'-f) with re- 
spect to (3 and 7, and inserting 7 = a and /3 = 1 to the 
obtained correlations to derive the final results. Once 
again, this can be done in a simple way for small lad- 
ders. Correlation functions (o'fj_;^(Tfj) and {<^2i-i'^2i+i) 
for spin ladder (|2.ip at a = 1 (quantum compass ladder) 
are shown in Fig. [7] for increasing temperature T. Other 
nearest neighbor correlations vanish at a = 1 for trivial 
reasons. 

Fig. [7] shows the qualitative difference between corre- 
lation functions of spin ladder (|2.ip and those of periodic 
QIM chain (|3.ip of length N, that appears in the ground 
subspaces r,; = r^+i. When all the subspaces are con- 
sidered, thermal fluctuations gradually destroy the spin 
order along the legs and the {o^i-iO^i+i) correlations 
weaken. On the contrary, the ((y^i-i'^^il correlations on 
the rungs are robust in the entire range of physically 
interesting temperatures < T < 2J, as the ZZ inter- 
actions destroying them are gradually suppressed with 
increasing T due to the increasing size of kinked areas. 

The above result is qualitatively different from the 
QIM results shown by dashed lines in Fig. [3 where ther- 
mal fluctuations initially increase intersite correlations of 
z-th spin components along the ladder legs and reduce 
the influence of the transverse field acting on rf pseu- 
dospins due to spin interactions 2 Jcfi.icrfj on the rungs. 
In the latter case thermal fluctuation in certain interval of 
temperature can enhance local spin ZZ correlations along 
the ladder legs at the cost of disorder in the direction of 
external field. This is because pseudospin interaction in- 
volves operators, not rf ones. Remarkably, in the full 
space, see solid lines in Fig. [71 the spin correlations are 
initially the same (at low T) as those for the QIM, but 
this changes when temperature T ~ 0.3J is reached and 
the two curves cross — then the rung correlations start 
to dominate. The crossing is caused by the growth of the 
kinked areas, as shown in Fig. [51 which are free of quan- 
tum fluctuations and therefore favor rung correlations of 
x-th spin components. 

Another interesting information on excitations in the 
quantum compass ladder is the evolution of the average 
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FIG. 7: (Color online) Nearest neighbor correlation func- 
tions, (crfi_ia|i} on the rungs (diamonds and red lines) and 
(crfi-icrfi+i) along the ladder legs (circles and blue lines), cal- 
culated for the compass ladder (a — 1) of 2N — 8 spins for in- 
creasing temperature T, taking into account: (i) all subspaces 
(solid lines) for increasing temperature T, and (ii) only the 
subspace which contains the ground state Vi = —1 (dashed 
lines). 

chain configuration with increasing temperature. As we 
know from Sec. IIV Al every subspace can be character- 
ized by the lengths of chains that appear in its label. 
Chain configurations can in turn be characterized by: (i) 
the number of chains which are separated by kinks to, 
and (ii) the total size of kinked areas K. Thermodynamic 
averages of both quantities, (m) and (K), can be easily 
determined at a = 1 even for a relatively large system 
using the combinatorial factor (|4.1ip calculated 

in Sec. IIVBI In the limit of T ^ oo one has: 

= ' ^'-'^ 

j:,L,}Fii{L^}M{L.}] 

where TO[{Li}] is the mmrber of {Li} in the chain config- 
uration Li, L2, • • • , L,n- 

In Fig. [5] we show the average quantities (m) and (K) 
for ladders of 2iV = 8 (left) and 2N = 104 spins (right). 
In both cases the average number of chains (m) starts 
from 1 and the average size of the kinked areas (K) starts 
from 0, corresponding to a single chain without kinks in 
the ground state at T = 0. The number of chains (m) 
grows to a broad maximum in the intermediate tempera- 
ture range and decreases asymptotically to a finite value. 
This non-monotonic behavior is due to the fact that the 
states with the highest energy, which become accessible 
when T — > 00, do not belong to the subspaces with large 
number of chains. The mean value of kinks (K) follows 
(to) but increases monotonically in the entire range of 
T, and for finite T one finds that (K) < (to). By look- 
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T/J 

FIG. 8: (Color online) Average size of the kinked areas {K) 
(|5.3p (circles and red lines) and the average number of chains 
(m) (|5.2|) (diamonds and blue lines) for the quantum com- 
pass ladder (|2.H) (at a = 1) consisting of: (a) 2N = 8, and 
(b) 2A*' — 104 spins. The mean size of kinked areas {K) in- 
creases monotonically with increasing temperature T to the 
asymptotic value N/4:, see Eq. (|5.4|l . The average number 
of chains (m) behaves differently, growing quickly to a maxi- 
mal value at intermediate T and then decreasing slowly when 
T ^ 00 towards: (a) 1.125, and (b) 13 + 12 x 10~^^ 

ing at the current results one may deduce that in case of 
T — > 00 and for large TV 3> 1 both quantities approach 

N 

(m)^ = (K)^ = -. (5.4) 

This is an interesting combinatorial feature of the chain 
configurations which is not obvious when we look at the 
explicit form of the combinatorial factor (|4.1ip . 
Note that Eq. (|5.4|) gives an integer due to our choice of 
system sizes 2N considered here, being multiplicities of 
8, i.e., TV is a multiplicity of 4. 



B. Spectrum of a large system 

The combinatorial factor Fi[{Li}] given by Eq. (|4.1ip 
enables us to calculate the partition function Z{1) (|4.4[) 
for a large system when a = 1. As a representative 
example we consider a ladder consisting of 2TV = 104 
spins. Even though we can reduce Hamiltonian (j2.ip to 
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a diagonal form when 2N = 104, as shown in previous 
paragraphs, it is still impossible to generate the full en- 
ergy spectrum for practical reasons — simply because the 
number of eigenstates is too large. Instead, we can obtain 
the density of states in case of a = 1 using the known 
form of the partition function (|4.4p and of the combina- 
torial factor (|4.1ip . Partition function for imaginary 1 /T 
can be written as 

Z{ix) = Yl = / dEe-'^'^piE) , (5.5) 

p=0 J Eo-e 

where 

4"-l 

p{E) =Y,5{E- Ep) , (5.6) 

and where sum is over all cigenenergies Ep of the ladder. 
Parameter Eq is the energy of the ground state. Small 
and positive e is introduced to formally include ±Eq into 
integration interval. Here we used the fact that ladder's 
spectrum is symmetric around zero at the compass point 
a = 1 (see Fig. [2])- Function p{E) can be easily recog- 
nized as the density of states. 

Using X — 27:n/w in Eq. (|5.5p . with w — 2(|£'o| + e) 
standing for the length of the integration interval and n 
being integer, we easily recover the density of states p{E) 
(|5.6p in a form of the Fourier cosine expansion 

P(E) = ^ E ^ (2-3 (2-£^) + ^^(0) ' (5-7) 

n— 1 

with amplitudes given by the partition function Z{ix). In 
practice we cannot execute the sum above up to infinity. 
Therefore, it is convenient to define Pc{E) which is given 
by the same Eq. (|5.7p as p{E) but where the sum has 
a cutoff for n = c. The heights of peaks in Pc{E) are 
expected to grow in an unlimited way with increasing 
value of c, so it is convenient to define the normalized 
density of states N{E) as 

N{E) = Pc{E)/pM . (5.8) 

The results for the compass ladder {a = 1) of 2N =104 
spins are shown in Fig. O These are relative density 
of states N{E) for cutoff c = 600 and Fourier coeffi- 
cients Z (2i7r^) for two intervals of n. Results obtained 
for lower cutoffs show that the overall gaussian shape of 
N{E), shown in Fig. [H^a), docs not change visibly if 
only c > 8. This allows us to conclude that the spec- 
trum of the compass ladder becomes continuous when 
the size of the systems increases which is not the case 
for the Ising ladder {a = 0). Higher values of n are in- 
vestigated to search for more subtle effects than gaussian 
behavior of N{E). These are found by looking at the 
amplitudes Z (2i7r^) in high n regime [Fig. [DJc)], as the 
low n regime [Fig. [D^b)] encodes only the gaussian char- 
acteristic of the spectrum. One finds three sharp max- 
ima of the amplitudes for n = 208, 330, 533 out of which 
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FIG. 9: (Color online) Relative density of states N{E) (a) 
as a function of energy ratio E/Eq and Fourier coefficients 
Z (2i7r;^) as functions of n for < ?i < 12 (b) and for 
13 < n < 600 (c) calculated for the ladder of 2N = 104 
spins. Relative density of states reminds a gaussian centered 
in zero with the width being roughly 0.15 of the spectrum 
width w. This follows from the gaussian behavior of Z (i27r^) 
coefficients for small n (panel (b)). Plot (c) reveals peaks 
in Z (i27r^) for n = 208,330,533, three order of magnitude 
weaker than for n = 0, corresponding with periodic conden- 
sations of the energy levels every AE — 1.28, 0.81, 0.50J (es- 
pecially every 0.81J). 



the one with n ~ 330 is about five times more intense 
than the rest, but it is still lO'^ times weaker than the 
peak in n = 0. These values of n correspond with some 
periodic condensations of the energy levels with periods 
AE = 1.28, 0.81, 0.50 J respectively which are visible in 
N{E) only in vicinity of £^ = ±£'o. 



VI. HEAT CAPACITY 
A. From Ising to compass model 

In this Section we analyze heat capacity to identify 
characteristic excitation energies in the compass ladder. 
We begin with complete results for the ladder consisting 
of 2N = 8 spins shown in Fig. [1] where all chain con- 
figurations can be written explicitly. Using Eq. (|4.4p for 
the partition function, one can next calculate all thermo- 
dynamic functions including average internal energy and 
the heat capacity. 

Results for the heat capacity Cy for different values 
of a arc shown in Fig. [TOl These plots cover three 
characteristic intervals of a where the behavior of curves 
changes qualitatively by appearance or disappearance of 
certain maxima. The positions of these maxima corre- 
spond to possible excitation energy scales of the system 
that change at increasing a and their intensities reflect 
the number of possible excitations in a given energy in- 
terval. In case of a = [Fig. fTUT a)]. we see a single 
maximum at ~ 2.2J which corresponds to flipping spins 
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FIG. 10: (Color online) Heat capacities Cv for spin ladder 
(|2.ip of 27V — 8 spins, shown in Fig. [T] with parameter a equal 
to: (a) a = 0,0.49,0.69,0.85, (b) a = 0.87,0.90,0.94,0.97, 
and (c) a = 0.982, 0.988, 1. In panels (a) and (b) lines from 
right to the left (solid, dashed, long-dashed, dashed-dotted) 
correspond with growing a. In panels (c) and (d) the values 
of Cv for a — 0.982, 0.988, 1 are shown by dashed, solid and 
dashed-dotted lines, respectively. Panel (d) shows the low 
temperature data of panel (c) (for a > 0.98), with a well 
developed small peak at low temperature originating from 
the critical excitations between subspaces Vi = 1 and = — 1 
close to a = 1; it disappears at a = 1. 



spends with low-energy mode of QIM that disappears for 
certain a > 1. This is not the case for other subspaces 
where QIM chains are fragmented and kinked area are 
formed. In case of the 1111 subspacc the low-energy 
peak in Cy vanishes at a ~ 0.65 and the high-energy 
peak persists and moves to higher temperatures with the 
increase of a. The situation is similar for the 1111 sub- 
space but the peak disappears at a ~ 0.75 and in the 
classical subspace 1111 we have only one maximum for 
any a. One can deduce now that the general rule is that 
the separation of peaks in heat capacity is reduced pri- 
marily by the growth of kinked areas and secondarily by 
the fragmentation of chains. This separation of energy 
scales is also visible in Fig. [2] where the spectra in differ- 
ent subspaces arc shown; below certain a in all cases but 
(d); which is the classical subspace, the energy gap be- 
tween the ground state and first excited state is smaller 
than other energy gaps appearing in the subspace. 

The mixing of different subspaces in the partition func- 
tion makes the peaks in Cy overlap which can result in 
reducing their number. This happens in Fig. Iior b): 
for solid (a = 0.87) and dashed {a = 0.90) curve we 
have only one maximum. For higher or lower a the en- 
ergy scales remain separated which is due to fact that: 
{i) soft modes survive in most of subspaces for low a, 
and (m) for high a the high-energy modes become even 
tougher and do not overlap with soft modes still present 
in subspaces with small kinked areas. The last phe- 
nomenon characteristic for the ladder are excitations be- 
tween Ti = —1 and ri = 1 subspace in the vicinity of the 
QPT. This yields to the appearance of the new energy 
scale A(q;) = 47VJ(1 — a) at a ~ 0.987 which manifests 
itself as a small peak in heat capacity in low tempera- 
ture. This maximum vanishes at a = 1 , as shown in Fig. 



in an Ising spin ladder. Switching on the XX interactions 
and weakening the ZZ interactions on the rungs has two 
effects: (i) decreasing energy and intensities of the high- 
energy maximum, and (m) appearance of a low-energy 
mode in every subspace with QIM chains which mani- 
fests itself as a peak with low intensity at low tempera- 
ture T, see Fig. [TUf a). At a ~ 0.85 this mode overlaps 
with modes of higher energies and until a ~ 0.94 there 
is a single peak again with a shoulder at high values of 
T, shown in Fig. [TOFb). Then the excitation energies 
separate again and a broad peak appears for high T ac- 
companied by a distinct maximum at T ~ 0.4 J. 

In Fig. [To] we recognize the characteristic features for 
the QIM chains present in most of the subspaces which 
are influenced by the excitations mixing different sub- 
spaces. If we had only one subspace with = —1, i.e., 
the one containing the ground state, then we would have 
two maxima in Cy for all < a < I — one of low inten- 
sity in the regime of low temperature T, and another one 
in high T, broad and intense. The small maximum corre- 



B. Generic features at large N 

After understanding the heat capacity in a small sys- 
tem of iV = 8 spin (Sec. I VI Ap . we analyze a large system 
using the statistical analysis of Sec. VI. Obtaining combi- 
natorial factor Fa [{Li}, R{Li}] in case of a < 1 is difficult 
and likely even impossible in a general way without fix- 
ing TV. Hence we focus on the compass ladder (a = 1) 
For the compass ladder of 2iV = 104 spins considered in 
Sec. IV Bl one finds 2^^ invariant subspaces. Although 
the eigenvalues can be found in each subspace, it is not 
possible to sum up over all subspaces for practical rea- 
sons and a statistical analysis is necessary. Therefore, 
the knowledge of the combinatorial factor Fi[{Li}], see 
Eq. (|4.1ip . is crucial to calculate partition function Z{1) 
(|4.4p . Fortunately, knowing it we only need to consider 
different chain configurations which are not very numer- 
ous — there are only 140854 of them. This means that 
on average each energy spectrum of the effective Hamil- 
tonian repeats itself almost 32 x 10^ times throughout all 
subspaces. 
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FIG. 11: (Color online) Heat capacity Cv for the compass 
ladder of 2N = 104 spins (spin ladder (|2.1|l at a = 1, solid 
line) as a function of temperature T . The main difference 
with the case of 2N = 8 spins [see Fig. llOf c)] is a small 
maximum appearing at very low T ~ 0.02J, see inset. This 
peak originates from the low-energy modes in subspaces Vi = 
Vi+i which exist in sufficiently long chains described by the 
QIM. Dashed line show heat capacity of the Ising ladder {a — 
0) of the same size. 



The statistical analysis of the compass ladder consist- 
ing of 2N = 104 spins in terms of: (i) mean values of 
kinked areas (K) (|5.3p . and {ii) the number of chains 
(m) (|5.2p . was already presented in Fig. [HJb), while the 
energy spectrum was discussed in Sec. IV Bl Here we 
present the heat capacity C'v for the compass ladder of 
this size in Fig. 1111 At high temperature one finds a 
broad maximum centered at T ~ 2J which originates 
from dense excitation spectrum at the compass point 
(a = 1). cf. the spectrum of the compass ladder with 
2N = 8 spins shown in Fig. [51 We remark that the 
broad maximum of Fig. [TT] has some similarity to broad 
maxima found in the specific heat (heat capacity) of spin 
glassesi^ However, here the broad maximum in the heat 
capacity does not originate from disorder but solely indi- 
cates frustration, similar as in some other models with 
frustrated spin interactions ji^^i We emphasize that the 
present results could be obtained only by developing a 
combinatorial analysis of a very large number of possible 
configurations of spin ladder, and due to the vanishing 
constant Cf{a = 1) = (|2.6p in the energy spectrum 
for the compass ladder. Unfortunately, the present prob- 
lem is rather complex due to the quantum nature of spin 
interactions, but in case of the binomial 2D Ising spin 
glass an exact algorithm to compute the degeneracies of 



the excited states could be developed recentlyi^ 

The heat capacity Cy of Fig. [Tl] at low temperature 
is qualitatively similar to the one obtained for 2N = 8 
spins, see Fig. llOf c). but the steep maximum at low 
T is here moved to lower temperature T ~ 0.2J. We 
also identified an additional (third) peak in the regime 
of rather low temperature T ~ 0.02J (shown in the in- 
set). This maximum originates from the QIM (j3.1|) where 
the energies of the ground state and of the first excited 
state approach each other for increasing N, if only a < 1- 
Thus, this lowest peak in the heat capacity obtained for 
the compass ladder of 27V = 104 spins has to be consid- 
ered as a finite size effect — for increasing system size 
it is shifted to to still lower temperature T, and would 
disappear in the thermodynamic limit TV — > oo, in agree- 
ment with the qualitative change of low energy spectrum 
of the QIM. 



VII. SUMMARY AND CONCLUSIONS 

We have investigated an intriguing case of increasing 
frustration in a spin ladder l|2.ip which interpolates be- 
tween the (classical) Ising ladder and the frustrated com- 
pass ladder when the parameter a increases from a ~ 
to a = 1. The ground state of the ladder was solved 
exactly in the entire parameter range by mapping to the 
QIM, and we verified that frustrated interactions on a 
spin ladder generate a QPT at a = 1, when conflicting 
interactions ZZ along the ladder legs compete with 2XX 
ones along the rungs. At this point the spin correlations 
on the rungs {(^2i-i'^2i) = ~^ collapses to zero and the 
ground state becomes disordered. We have shown that 
the ground state of a finite ladder has then degeneracy 2, 
while the analysis of the energy spectra for increasing size 
suggests that the degeneracy increases to 4 in the ther- 
modynamic limit. We note that this result agrees with 
degeneracy 2x2^ found for the 2D compass model,— 
where L is a linear dimension (the number of bonds along 
one lattice direction) of an L x L cluster in the 2D sys- 
tem. In our case of a 2 x iV ladder, L = 1 for ladder 
rungs, so indeed the degeneracy is 2 x 2 = 4. 

The present method of solving the energy spectrum in 
different subspaces separately elucidates the origin of the 
QPT found in the present spin ladder (|2.ip at the point 
a = 1 , corresponding to the frustrated interactions in the 
compass ladder. We argue that this approach could help 
to find exact solutions in a class of quasi- ID models with 
frustrated spin interactions, but in some cases only the 
ground state and not the full spectrum can be rigorously 
determined. For instance, this applies to a spin ladder 
with frustrated spin interactions between different triplet 
components on the rungsj^ where a different type of a 
QPT was found recently. 

By performing a statistical analysis of different possible 
configurations of spin ladder (j2.ip with periodic bound- 
ary conditions we derived a partition function Z{a) for 
a mesoscopic system of 104 spins. The calculation in- 
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volvcs the classification of ladder subspaccs into classes 
of chain configurations {Li} equivalent by symmetry op- 
erations and the determination of the combinatorial fac- 
tor Fa[{Li}, R^j^. j]. We have shown that this factor can 
be easily determined at the compass point {a = 1), so 
the heat capacity of such a mesoscopic compass ladder 
could be found. 

Summarizing, we demonstrated that spin ladder stud- 
ied in this paper exhibits a QPT from a classical ordered 
to a quantum disordered ground state which occurs due 
to the level crossing, and is therefore of first order. It 
leads to a discontinuous change of spin correlations on 
the rungs when the interactions along the ladder legs 
and on the rungs become frustrated. Fortunately, the 
subspaces which are relevant for the QPT in the com- 
pass ladder considered here can be analyzed rigorously, 
which gives both the energy spectra and spin correla- 
tion functions by mapping the ladder on the quantum 
Ising model. The partition function derived in this work 
made it possible to identify the characteristic scales of 
excitation energies by evaluating the heat capacity for a 
mesoscopic system. 

Note added in proof. After this paper was accepted, 
we learned about a powerful algebraic method to ana- 
lyze exactly solvable spin Hamiltonians.— The present 
quantum compass laddeed could be also analyzed using 
this approach. 
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APPENDIX: EVALUATION OF THE ENERGY 
ORIGIN Cp(q) in invariant SUBSPACES 

We need to express X^iLi ''i' which appears in Cp{a), 
see Eq. (|2.6p . in terms of chain configurations {Li}. This 
task may be accomplished by the following construction. 



Let us imagine certain string of r^'s written in terms of 
chains {Li} and kink areas {^i}: 

Ai{Li)A2{L2)A3{L3) ■ ■ ■ Ak{Lk) . (A.l) 

First, we want to calculate the sum of r,;'s included in 
chains. We choose any from the chain Li and fix its 
sign as Hn. Now this chain gives rij^Li contribution to 
the total sum of r^'s. To get to the second chain we have 
to pass through the first kink area Ai. If the number 
of kinks in Ai is even, then the next chain will give the 
contribution rinL2, and if not, then it will give the oppo- 
site number. Therefore, after passing through the whole 
system we will get the term 

ri„(Li + P2L2 + P2P3L3 -I- ... -I- P2P3 ■ ■ -PkLk) , (A. 2) 
where pi = (—1)^% and Ki is a number of kinks in 
kink area i. It is clear that the parameters {pi} sat- 
isfy HiLiPi = 1- Now we need to calculate the sum of 
r-i's placed in kink areas. The sign of the first chain is 
already chosen as rin so we pass to A2. For even number 
of kinks in A2 the contribution is zero. If the number is 
odd, then we get the sum equal — rin. Passing to the next 
kink area we follow the same rules but wc have to change 
Tin into 752^11 • The total contribution from the kink areas 
is then equal to 

jl+Pl A I+kI q^ 

-Pir-mi — \- 2_^piP2. ■ - Pr-i — - — > • (A.3) 

Using the results given in Eqs. (1X211 and (jXSl) we obtain 
finally 

N r k ^ 

^n = rinlLi-l + ^P2P3...Pi{L,-l) \ . (A.4) 

i=l I i=2 J 

Thanks to this result, wc can write the energy given by 
Eq. (|4.3p in terms of new variables {Li,pi} instead of 
{r,;} which arc definitely more natural for the present 
problem. 
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